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Abstract 

Shape preservation behavior of a spline consists of criterial conditions for preserving 
convexity, inflection, collinearity, torsion and coplanarity shapes of data polgonal arc. We 
present our results which acts as an improvement in the definitions of and provide geomet- 
rical insight into each of the above shape preservation criteria. We also investigate the effect 
of various results from the literature on various shape preservation criteria. These results 
have not been earlier refered in the context of shape preservation behaviour of splines. We 
point out that each curve segment need to satisfy more than one shape preservation criteria. 
We investigate the conflict between different shape preservation criteria l)on each curve 
segment and 2)of adjacent curve segments. We derive simplified formula for shape preser- 
vation criteria for cubic curve segments. We study the shape preservation behavior of cubic 
Catmull-Rom splines and see that, though being very simple spline curve, it indeed satisfy 
all the shape preservation criteria. 



1 Introduction 

Designers in industries need to create splines which can interpolate the data points in such a 
way that they preserve the shape of polygonal arc formed by data points. Among the properties 
that the spline curves need to satisfy following properties are of common interest to almost all 
the designers: 

• Smoothness 
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• Preservation of shape of the data polygon 

• Each curve segment to be a low order polynomial curve. 

We first illustrate the shape preservation behaviour of a spline interpolating planar data 
points with the help of figures [fl [2] and [3] 




Figure 1: Data points to be interpolated 





Figure 3: Data points with shape preserving interpolating spline 

• Inflection depicted by data points x 3 x 4 x 5 and x 6 is preserved by suitable inflection of 
the curve segment between x 4 and x 5 

• Convexity depicted by data points x 4 x 5 x 6 and x 7 is preserved by the convex shape of 
the curve segment between x 5 and x 6 

• Collinearity depicted by data points x Xi and x 2 is preserved by collinearity of the curve 
segment between x and x 2 

an so on. It can also be observed that shape preserving behaviour of a spline makes it more 
close to mimicing free-hand curve drawing. 

However, modelling of shape preservation behaviour of interpolating splines in R 3 is rela- 
tively difficult. Typically shape preservation criteria that have been studied for the generation 
of interpolating splines consists of conditions for the preservation of l)Convexity 2)Inflection 
3)Collinearity 4)Torsion and 5)Coplanarity shapes of data polygonal arc (formed by line join- 
ing the consecutive points of ordered set of data points). In this paper we present literature 
survey and also our analysis for each of the above criteria. In sections [3TI and l3l we present our 
analysis for convexity preservation criteria. In literature, basically, two definitions for convexity 
preservation criteria are followed. One which is followed in [8, Kaklis and Karavelas, 1997], 
lfT0l[T2l Costantini et. al.] etc. and another which is followed in lfT3l Kong and Ong, 2002], [[71 
Goodman and Ong, 1997] etc. The latter definition includes the conditions of previous defini- 
tion and hence we investigate further on the later definition. We find that in the later definition, 

3 



for a large set of interpolating spline curve (containing set of rational spline, excluding straight 
lines, curves as small subset) one of the two conditions is redundant (that is, the condition is 
actually taken care by other condition of the definition) for all the points of the curve except for 
set of points on the curve of measure zero (the set is finite for rational curves excluding straight 
lines). Also at these points the error (if at all occurs) is very negligible and it is observed that 
in all almost all the algorithm staright line segments in the spline curves are considered for 
collinearity preservation criteria. We use our lemma I3T31 and the characterization of convexity 
of planar curve presented in [[61 Liu and Traas, 1997] for our analysis of convexity criteria of in- 
terpolating splines. In j6l Liu and Traas, 1997] the characterization of convexity of planar curve 
on R 2 , that is, XY plane have been derived. We use this characterization to get the characteri- 
zation of planar curve on any plane in R 3 . Then the modified characterization has been used to 
improve the definition of convexity preservation criteria of splines. We further state simplified 
characterization of convexity preservation criteria for cubic splines in terms of control polygon 
of individual Bezier segments. 

In sections |4] to [6] we present the analysis for inflection preservation criteria of splines. We 
refer two papers [5, Goodman, 1991] and Goodman and Ong, 1997] for our analysis for 
inflection criteria. In [5, Goodman, 1991] author has stated definitions and results for inflection 
counts for planar and space curves and polygonal arcs. In |[7l Goodman and Ong, 1997] au- 
thors have defined the inflection criteria of splines and have constructed a spline satisfying the 
criteria. However, in [[7J Goodman and Ong, 1997], authors haven't indicated any connection 
with the analysis of J5] Goodman, 1991]. Also the analysis in [5, Goodman, 1991] is not used 
in |[7l Goodman and Ong, 1997] for the analysis of inflection criteria of splines. We observe in 
section [6] that our lemma 137131 acts as a connection between the analysis in sections 0] and \5\ and 
conditions stated in definition of convexity criteria of splines. 

In section|4]we state definitions and results from [51 Goodman, 1991] for inflection count for 
planar curves and polygonal arcs. The relation between the inflection counts of planar B-spline 
and Bezier curves and inflection counts of their control polygons are stated. In section [5] we 
state definitions and results for inflection counts for space curves and polygonal arcs from [[51 
Goodman, 1991]. Definitions in section \5\ uses the definitions in sectional In section[6]we state 
the definition and analysis for the convexity criteria for splines. In section[7]we state the results 
from [5, Goodman, 1991] which states (via lemma [37131) the conditions under which a spline 
curve does not satisfy inflection criteria or convexity criteria. 
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In sections [8] [9] and [TO] we analyze and give improved conditons for collinearity , torsion 
and coplanarity preservation criteria respectively for splines. In section [TTJ and [12] we analyze 
and give condtions on a spline curve to resolve conflict between different shape preservation 
conditions on a curve segment or a pair of adjacent curve segments. 

Almost every geometric modeler necessarily uses cubic splines for generating 3D models 
of products. Cubic splines are computationally most viable solution for various applications 
requiring complicated geometric operations. They are also the splines of least degree that can 
exhibit torsion. Thus we see that it is necessary that cubic splines should preserve shape of the 
data points they interpolate. In section [13] we describe our analysis and results for shape pre- 
serving criteria for cubic splines. The conditions for shape preserving criteria for cubic splines 
derived in this section are expressed in terms data points and slopes at data points. During our 
analysis we obtain simplified formula for discrete shape measures and a new property for Bezier 
curves. 

In subsection 113.11 we state aome additional notations required in section [13] In subsec- 
tions [1321 [1331 [EH [[33] and \TJM we derive the conditions for convexity, inflection, torsion, 
collinearity and coplanarity preservation criteria respectively, for the cubic curve segment in 
terms of the two data points at their ends and slope vectors at them. In subsection 113.41 we 
simplify the expression for torsion of cubic Bezier curves. In subsection [133] we derive an ex- 
pression for sine of the angle between a point vector on a Bezier curve with a given vector and 
state it in theorem 113.171 and we use it to get simplified formula for collinearity preservation 
criteria for cubic curve segments. Further we use this analysis in subsection 113.61 to get sim- 
plified formula for coplanrity preservation criteria for cubic curve segments. In section [T4l we 
investigate shape preservation behavior of cubic Catmull-Rom splines. Finally in section [T5l we 
state our conclusions about the analysis in this paper. 



2 Notations and prelimnaries 

Let Xj G R 3 , i = 0, n be n + 1 data points and V be the polyline or polygonal arc formed 

by joining the points with each side being Li = Xj — Xj_ l5 i = 1, n. Let iVj = Lj_i x Lj. In 

_ Ni 
discrete differential geometry Hi Sauer, 1970] discrete binormal is defined as rjyl- F° r a curve 

7 (t) = [x(t),y(t),z(t)],te [0,1] ini? 3 let^(t) = 7 '(t) x 7 "(t). 
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3 Convexity preservation criteria for interpolating splines 

In E Karvelas and Kaklis, 2000], © Kaklis and Karavelas, 1997] El Costantini, Goodman, 
Manni, 2000] EH Costantini, Cravero, Manni, 2002], 031 Manni, Pelosi, 2004] we have the 
following definition 

Definition 3.1 Convexity preservation criteria for a curve j(t) interpolating data points con- 
sists of following condition: 

1. ifNi-i ■ Ni > 0, then uj(t) ■ N m > 0, t G U], m = i-l,i. 

The above definition 13.11 takes into account the ability of a spline to appear as convex curve 
along only two viewpoints AT i _ 1 and iV,. 




Figure 4: Data point with N^i ■ N { > requiring convexity preservation by i curve segment 

The following definition is from [7, Goodman and Ong, 1997], Q~3l Kong and Ong, 2002]. 

Definition 3.2 /[7| Goodman and Ong, 1997] Convexity preservation criteria for a curve j(t) 
interpolating data points consists of following condition: 

1. IfNi-x ■ N l > 0, then for all N = XN^ + fiN it where X,p > 0, (X, p,)^(0,0), the 
projection P N ±j(t), t G \pi-i, U], is globally convex and 

2. u(t) ■ N > 0. 

Now we state a theorem from [JVl Goodman and Ong, 1997] which makes the convexity 
condition of projection curves simpler. 

Theorem 3.3 /:7 Goodman and Ong, 1997] Suppose that R is a curve in R 3 and V\ and V% are 
vectors in R 3 so that P v ± R and P v ± R are convex with the same orientation with respect to V\ 
and V2 respectively. Then Py±R is convex with the same orientation with respect to V where 
V = \V X + fiV 2 for any A, fi > 0, (X, fi)^(0, 0). 
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Thus from the theorem (13.31 ) we can observe that condition requiring P N ± 7, (t) and P N ± 7; (t) 
to be convex is equivalent to the condition requiring curves P N ±ji(t) for N = AiVj_i + /xiVj for 
any A, p, > 0, (A, p)^(0, 0), to be convex. 

3.1 Convexity of planar curves 

3.1.1 Convexity of planar curves on XY plane 

In [0 Liu, Traas, 1997] authors have defined local and global convexity of a planar curve as 
follows. The author distiguishes convexity and concavity of planar curves in terms of the orien- 
tation we assign to the curve. Consider a curve 7(t) = (x(t),y(t)), t G [0, 1] in R 2 . An oriented 
planar curve is an ordered set in R 2 , given by 7(2) = [x(t),y(t)], t G [0, 1] with direction from 
t = to t = 1. A global supporting line of an oriented curve 7(2) at a point 7(^0) is an oriented 
line, L, having consistent direction with 7(2) in t > and satisfying (a) 7(^0 ) is a point of L; (b) 
the entire curve 7(i), i G [0, 1], lies in one closed half-plane with respect to L. 

A global supporting line of an oriented curve j(t) at a point 7(^0) is an oriented line, L, 
having consistent direction with 7(2) in to, and satisfying 

1. 7(io) is a point of L; 

2. The entire curve 7(t), £ G [0, 1], lies in one closed half-plane with respect to L. 

A local supporting line of an oriented curve j(t) at a point 7(^0) is an oriented line, L, 
having consistent direction with 7(2) in to, and satisfying 

1. 7(io) is a point of L; 

2. A local neighborhood 7(t), t G [ti, t 2 ], of 7(io)> li es in one closed half-plane with respect 
to L, where ti and t 2 satisfy 

< h < t < t 2 < 1 or = t x = to < t 2 < 1 or < t x < t = t 2 = 1. (3.1) 

Definition 3.4 Liu, Traas, 1997] j(t), t G [0, 1], is a globally convex curve if it satisfy: 

1. There is at-least one global supporting line at every point ofj(t); 

2. The entire curve lies in the right closed half-plane with the supporting line as its left 
boundary. 
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Figure 5: Supporting and non- supporting lines [Liu and Traas '97] 





(a) closed (b) open 

Figure 6: Globally convex curve [Liu and Traas '97] 

Definition 3.5 /|6] Liu, Traas, 1997] j(t), t e [0, 1], in definition di.4l) is a locally convex curve 
if in (1) the global supporting line is replaced by a local one and (2) is true only for a related 
local neighborhood. 

A globally and locally concave curve has the same definition as for globally and locally 
convex curve respectively with left and right interchanged. 

We observe that the above definitions [J hold for any planar curve in R 3 . However, the two 
main theorems of (6] Liu and Traas, 1997] that we state below need some modifications. 



'There are other definitions of a convex curve. For example, in J4] Farin] it is defined as a part of boundary of 
a convex set. In Liu and Traas, 1997] is has been proved, using Hahn-Banach theorem, that this definition is 
included in definition ( 13.4b . 
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(a) 




Figure 7: (a) Locally convex curve [Liu and Traas '97] (b) Condition 13 .61 for global convexity 
[Liu and Traas '97] 



Theorem 3.6 ^ Liu and Traas, 1997] j(t), t G [0, 1], is locally convex if and only if 

j'(t) x i\t) <o,te [o, l] 

where 7 (t) is C 2 -continuous, j'(t), j"(t) are first and second derivatives ofj(t), and 



i{t) x i\t) 



x'(t) x"(t) 
y'{t) y"(t) 



x'(t)y"(t)-x"(t)y'(t). 



(3.2) 



(3.3) 



Theorem 3.7 jj^ Liu and Traas, 1997] A curve ■jit) satisfying the condition 



T (t)^ 7 (0)/orte(0,l) 



(3.4) 



is globally convex if and only if 



7 '(t)x 7 "(t) < 0, 

( 7 (t)-7(o))xy(t) < at g [o,i], 
y(o) x ( 7 (t) - 7 (o)) < o,te[o,i]. 



(3.5) 
(3.6) 
(3.7) 



Remark 3.8 For a curve j(t) G R 2 , j'{t) x j"(t) is its curvature. Therefore, in \3.2\ and \3.5\ 
the equality holds att — t\ if and only ifj(t) behaves locally as straight line (turns with angle 
0°) at t = t\. One can replace \3~2\ and [5731 by j'{t) x 7 "(t) < by requiring strict convexity 
ofjif). But we did not find this to be significant requirement as set of points for which equality 



holds is of measure zero for almost all curves (except straight lines) used for interpolation (and 
therefore their local behaviour as straight line have negligible effect). The conditions \3. 6\ and 
\3.7\ ensures that the curve doesn 't intersect itself and equality in these conditions may cause the 
curve to become straight line in some cases. 




(a) (b) 
Figure 8: Non-convex curve (a) without [331 (b) without [377] [Liu and Traas '97] 



3.1.2 Convexity of planar curves in R 3 

In this subsection we extend the characterization of globally convexity of curve in R 2 to the 
characterization of global convexity of planar curve j(t) = [x(t), y(t), z(t)\, t 6 [0, 1] in R 3 . In 
the remaining part of this chapter, except in section 0] we use the operator x to denote the cross 
product of two vectors, that is, 



Ax B 



A X Ay 

B,,. B„, 



k 

A z 
B~ 



(3.8) 



where A = [A X) A y , A z ], B = [B x , B yi B z ] e R 3 , i, j and k are unit vectors along x, y and z 
axis respectively. We observe that [6, Liu and Traas, 1997] the operator x (in R 2 ) is mainly used 
to understand the direction in which the curve bends with respect to the orientation induced by 
z— axis. We now explain our characterization for global convexity of planar curve in R 3 which 
solves our above purpose as follows. 

It has been observed in [[6], Liu and Traas, 1997], the convexity and concavity depends on its 
orientation, that is, direction in which is traversed. By inverting the orientation, a convex curve 
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turns into a concave curve and vice versa. Definition based on direction helps to distinguish 
between two curves as convex and concave curves in most of the practical situations. For 
example, in the case of two curves y(t) and Q(t) (t e [0, 1]) such that 7(1) and Q(l) lie on the 
right side of 7(0) and Q(0) respectively. 




Figure 9: Convex curves with different orientations [Liu and Traas '97] 

But the direction of orientation itself gets inverted if the normal to the plane is taken to be the 
normal with its direction opposite to the given normal. (Note that in (|3.3I) " x " can be interpreted 
as to denote the dot product of the cross product between two vectors in x — y plane with unit 
vector along z-axis (which is along normal dierction to the x — y plane)). In our case we need 
the convexity of the curve according to orientation induced by a specified normal vector as it 
requires that the spline curve to be convex in the same direction as the (data) polygonal arc. The 
normal vector is specified as the normal to the plane containing a pair of adjacent line segments 
in the data polygonal arc. (Therefore we do not consider the case of concavity.) 

Thus we have the following definition for local and global convexity of a planar curve lying 
on a plane II in R 3 with respect to orientation of normal vector iV of the plane IT. 

Definition 3.9 A C 2 -continuous planar curve j(t), t £ [0, 1], is locally convex if and only if 
(y'(t)xy"(t))-N>0,te[0,l]. 
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Definition 3.10 A C 2 -continuous planar curve j(t), t G [0, 1], satisfying the condition 

7 (*)^7(0)/or*e(0,l) (3.9) 
is globally convex if and only if it satisfies following conditions. 

1. ( 7 '(t) x Y(t)) ■ N > 

2. (( 7 (t)-7(0))x 7 '(t))-iV>atG[0,l] 

3. (y(o)x (7(0-7(0))) -iv>ate [0,1] 

Following remarks about defintions 13.91 and [3 . 1 01 are very significant. 

Remark 3.11 Since 7'(t), 7"(f) and ("f(t) — 7(0)), t e [0, 1] are parallel to the plane IT. 
Therefore sign of dot product of each cross product with N actually represents the bending of 
the curve according to the orientation induced by N. 

Remark 3.12 Remark \3~M about equality holds for the inequalities appearing in the definitions 
UMandUUhholds. 

3.2 Improvement in the condition of convexity preservation criteria 

Now using our lemma, stated below, we find that for almost all interpolating splines the condi- 
tion u(t) ■ N > in definition 13 .21 is implied by global convexity of P N ±( / y(t)) (and need not 
be stated separately) for almost all values of t. Thus we further simplify definition of convexity 
criteria. This lemma also helps to modify the definition for inflection criteria to get a simpler 
definition in section [6l 

Lemma 3.13 Let P N ±{~f(t)) be denoted as 7iv(t), u(t) = 7'(t) x j"(t) and u N (t) = "/ N (t) x 
7&(t). Then u N (t) ■ N = oj(t) ■ N. 

(x, y, z) • N + d 

Proof: A plane with normal vector N is given by — - — - = 0, d 6 R . We know that 



/N 7 (t) ■ N + d^ T m p ... .. , -y'(t) -N + d^ 
lN (t) = 7(f) + ,, A n|2 N - Therefore i N {t) = 7 '(f) + mr N > 



n/'(A . N + d 

7* (*) = 7"(f) + ,| An , 2 N. From the above we get 



Y(t)-N + d^ t f i(t)-N + d 
}N 

u N {t) ■ N = 



UN (t) = ( 7 '(f) x 7 "(f)) + (7'(f) x 1 N) + ( '"a N x /(f)). Thus 
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JJV,, 

Hence proved. I 
Now we first analyse the condition u>(t) ■ N > 0. 

• For a large class of curves, including rational curves (except for straight lines) as a small 
subset, u(t) ■ N = holds for a set of values of t G [0, 1] whose measure is zero. 

• uj(t) ■ N = at t = t\ if and only if |o>jv(t)| = that is, 7 at(£) behaves as straight line at 
t = t\. Such behavior at t = t x and not in its neighbbourhood, has negligible effect on 
shape of the projected curve (along a viewpoint). 

• For some cases u(t) ■ N = at t = t\ may also imply that u(t), which is binormal of 
the curve j(t), is perpendicular to N at t = t\. That is, is parallel to the osculating 
plane (plane on which curve lies locally) of 7 (t) at t — t\. This might not be good unless 
such a torsion is required. Also, this may occur to great extent even for the case where 
e > u(t) ■ N > 0, for sufficiently small values of e. However, possibility of such torsion, if 
undesired, can be controlled by other shape preservation criteria like torsion, coplanarity 
and collinearity preservation criteria. 

From the above analysis and lemma I3T31 we get following two observations. First, from 
definition 13 .91 (for local convexity of a curve), we see that the definition 13 . II requires the projec- 
tion of curve on the plane perpendicular to and Ni be only locally convex which one can 



observe from figure [8(a)] that the spline may not always serve the purpose of shape preservation. 

Second, according to the definition 13.101 for convexity of a planar curve in R 3 one of the 
condition that 7jv(t) need to satisfy, to be globally convex, is co>jv(t) • > 0. From our lemma 
13.131 global convexity of P N xj(t) implies tu(t) ■ N > 0. Thus we see that the condition 
u>(t) ■ N > is redundant in the definition 13 .21 

We now state our results which will further simplify the conditions in the definition 13 .21 as 
lemmas. The proofs of these lemmas are similar to that of lemma I3T31 

Lemma 3.14 With the notation same as that in lemma \3J3\ we have 

((7^(t) - 7iv40)) x 7 ^(t)) • N = (( 7 (t) - 7 (0)) x 7 '(t)) • N 

Lemma 3.15 With the notation same as that in lemma \3J3\ we have 

(7^(0)x( 7iV x(0-7Arx(0)))-Ar = ( 7 (0) x ( 7 (t) -7(0))) -AT 
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Using theorem [331 and lemmas [3 . 13U3. 14| and l3 .151 we have the following improved defini- 
tion for convexity criteria for splines. 

Definition 3.16 A spline curve j(t) interpolating data points satisfies convexity criteria if, for 
j = % - l,i 

1. u(t) ■ Nj > 0, 

2. (( 7 (f) - 7 (0)) x Y(0) • Nj > 0, 

3. (Y(0)x( 7 (0-7(0)))-iV i >a 
t E [U-i, U], whenever iVj_i • Ni > 0. 

4 Inflection of a planar curves and polygonal arcs 

In [[5l Goodman, 1991] authors have given definitions and conditions for the existence of inflec- 
tions in a planar curve as follows. In this section for A = (A%, A 2 ), B = (Bi, B 2 ) in R 2 we 
write 

Ax B = A 1 B 2 -A 2 B 1 (4.1) 

For any sequence a = (ai, . . . , a n ) in R n , we define S(a) = S(a\, . . . , a n ) to be the number of 
strict sign changes in the sequence. 

Definition 4.1 /O Goodman, 1991 ] We say a polygonal arc PqP\ . . . P n for points Pq, P\, . . ., 
P n in R 2 is regular if the following hold 

1. It turns through a total angle of magnitude at most tt, that is, for some V in R 2 , V ■ (Pi — 
Pi-i) >0,i = l,...,n. 

2. It does not turn through an angle of tt at any vertex, that is, for any < % < j < n with 
Pi..! + Pi = P i+1 . . . = Pj + P J+l , Pi - Pi^ ^ \{P j+1 - Pj) for any A < 0. 

Definition 4.2 /[5] Goodman, 1991 ] For a regular polygonal arc PqPi • • • P n in R 2 , with the 
condition Pi-\ ^ Pi, i = 1, . . . , n denoted by F we define inflection count as 

i(r) = ^....K-i) (4.2) 

where 

Vi = (P i -P i _ 1 )x(P i+1 -P i )i = l,... 1 n-l. (4.3) 
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For any function / : (a,b) — > R we define S(f) to be the number of strict sign changes in 
f(t), a < t < b, that is S(f) = supS(f(ti), . . . , f(t n )), where the supremum is taken over all 
sequences a < t% < . . . < t n = b, for all n. For a curve j(t) E R 2 , t e [a, b] which is constant 
for t in [a, 0\ C (a, 6) but not on any larger interval, we define 

f Uu(a~) xu'(a-)+u(P~) xu'(p-)}, if u(cr) = u(p + ), 
K{t) — I 2 (4.4) 
(«(a-)xu(^), ifw(a-) ^ w(/J+). 

Definition 4.3 Inflection count for a curve jit) 6 i? 2 , t G [a, 6], is defined as i(j) := S{K). 
It is actually the number of times the curve changes from turning in a clockwise direction to 
turning in an anti-clockwise direction, or vice-versa. 

With the above definitions regarding inflection count we state the following relation between 
inflection count of B-spline curve and its control polygon from [5, Goodman, 1991]. 

Theorem 4.4 £5] Goodman, 1991 ] Suppose 

m—l 

^t) = A i N i(t)>to < t < tm, (4-5) 

i=—n 

and r denotes the polygonal arc A_ n A_ n+ i ■ ■ ■ A m _i. If Ai_ n ■ ■ ■ Ai is regular for i = 0, . . . , m— 
1, then 

tii) < tin (4.6) 

where N{ denote B-spline basis function. 

We now state the relation from [5 , Goodman 1991] between inflection count of Bezier curve 
and its control polygon which follows as a corollary to the above theorem. 

Corollary 4.5 /[5] Goodman, 1991 ] If 

n 

jit) = Y,M?)t l i~L-t) n ~ l ,0<t<l, (4.7) 
and the polygonal arc T = A A 1 ■ ■ ■ A n is regular, then 

iii) < i(r). (4.8) 



Now we state a theorem from [5, Goodman, 1991] which says that the above result does not 
hold if the control polygon is not regular. 
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Theorem 4.6 £5] Goodman 1991 ] Let j(t) be cubic Bezier curve given by 

7(t) = A(l-t) 3 + Wt(l-tf + 3Ct 2 (l-t) + Dt 3 (4.9) 

where A, B, C, D E R 2 are control points and Y denote the polygonal arc ABCD. Suppose 
r turns through an angle of magnitude > tt and let P be the point of intersection of the line 
through A and B and the line through C and D. Then 




(4.10) 



Figure 10: Cubic Bezier curve with two inflection points whose control polygon has no inflec- 
tion 



5 Inflections of curves and polygonal arcs in R 3 

Let 7 : [a, b] — > R 3 be a curve in R 3 . We also resume the meaning of operator x as cross 
product of two vectors in R 3 as defined in subsection 4.1.3. 

For any w in Si = {v E R 3 : \v \ = 1} we shall denote by P w the orthogonal projection 
from R 3 onto the 2-dimensional subspace orthogonal to w, that is, P w x = x — (x ■ w)w. 

Definition 5.1 /T5] Goodman, 1991] Inflection count 1(7) of the (spatial) curve 7 to be the 
maximum number of inflections that be seen in 7 by observing from any direction, that is, 

1(7) = ess sup{i(P w , y) : w E Si}. (5.1) 
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Figure 11: Curvature plot of the curve in figure [TOl 

We suppose that 7 is continuous with piecewise C 1 unit tangent vector u(t) = j'{t)/\j'(t)\. 
As before suppose 7(2) is constant for t in [a, (3} C (a, b) but not on any larger interval and we 
define for a < t < (3, 

/ \{u(ar) x u'ipr) + u(p + x u\l3+)), \fu{a~) = u(f3 + ), 

(5.2) 

u(a~) x u(f3 + ), if u(a~) ^ u(f3 + ). 

Theorem 5.2 £5] Goodman, 1991 ] Suppose 7 : [a, b] — ► i? 3 ?5 continuous with piecewise C 1 
unit tangent vector u(t) = , 7 ,[!L 77ze/7 



A-(t) 



7(7) = sifp {^(w.K) : iu G 5i} (5.3) 

Corollary 5.3 If the curve 7, a* in theorem 0.2D . //e* in a plane with a normal n, then 

J( T ) = S(n-K). (5.4) 

Theorem (15.21) . also implies the following definition for polygonal arc. 

Definition 5.4 /|5] Goodman, 1991] For a polygonal arc PqP\ ■ ■ ■ P n and Pi-\ ^ Pi, % — 
1 , . . . , n, denoted by V then its inflection count is 



I(T) = sup {S(w ■ Vi, . . . , w ■ V n ) : w G Si}. 



(5.5) 



6 Inflection preservation criteria for interpolating splines 



Let Xj G i? 3 , i = 0, n be n+ 1 data points and X> be the polyline joining the points with each 
side being Lj = — Xj_i, i — 1, n. Let iVj = Lj_i x Lj. In discrete differential geometry 

17 



[Sauer, 1970] discrete binormal is defined as Tph- F° r a curve P{t) = [%(t),y(t), z(t)], t G 
[0, 1] in R 3 let u(t) = P'(t) x P"(t). 

In |HSl Costantini, Goodman, Manni, 2000], OH Costantini Cravero Manni, 2002] lfT51 
Manni, Pelosi, 2004] we have the following definition 

Definition 6.1 470] Costantini, Goodman, Manni, 2000] A curve 7(i) interpolating data points 
satisfies inflection criteria if is satisfy the condition that, ifN^i ■ Ni < 0, (^(^) ■ N m )(Ni ■ 
N m ) > 0, l,m = i — and u(t) ■ Nj, has precisely one sign change in t G 
m — % — 1, %. 

The above definition requires the projection of curve on the plane perpendicular to iVj_! and 
Ni have only one inflection point. Thus it takes care about inflection preservation along two 
viewpoints only. 

In [|7l Goodman and Ong, 1997], HH Kong and Ong, 2002] we have the following definition 

Definition 6.2 /[7| Goodman and Ong, 1997] Inflection preservation criteria is defined by the 
condition that ifNi-i ■ Ni < 0, 

1. u(ti-x) ■ Ni-i > 0, uj(ti) ■ Ni > and 

2. for all N = XN^i + fiNi, where Xfi < 0, (X, /iJ^fO, 0), u(t) ■ N has precisely one sign 
change in tj]. 

One can see that the above definition takes care about inflection preservation along all the view- 
points between N^i and Ni along the plane containing the two normals. 

We now observe that our lemma [3T31 acts as the connection between the analysis of inflec- 
tion of curves and polygonal arcs and the definition l6.2l The two conditions basically states that 
the projection of the curve 7(t) on the plane with normal vector N, P N ±( / y(t)), t G 
should have only one inflection point. 

We now state the condition under which a curve segment of "f(t) with Bezier representation 
satisfies the inflection criteria. Let a Bezier curve j(t) be 

n 

7 (t) = ^^(™K(i-tr-\o<t<i, (6.D 

Ai G R 3 , the polygonal arc T = A Ax ■■■A n ,V l = {A, t - A^) x (A i+1 - A { ). 
If Ai, i — 0, n are such that 
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Figure 12: Data point with A^_i • Ni < requiring inflection preservation by i curve segment 

1. Ni ■ V\ < 0, Ni ■ V n -i > and P N ± have only one inflection 

2. N 2 ■ Vi > 0, N 2 ■ V n -i < and P N ± have only one inflection, 

3. (N 1 -V i )(N 2 -VJ<0,i = l,...,n-l. 

and we have two scalars A, // such that A/i < 0, (A, /i) 7^ (0, 0), iV = AA^j + /jN i+ i then we 
have Vi ■ N = XVi ■ iVj + /iV^ • A^ +1 and thus P^±T has only one inflection point. 

Thus we see that if the i th curve segment of interpolating spline ji(t) is a cubic curve and 
satisfies the first condition of the definition (16.21) then it also satisfies the second condition. 

7 Difficulties in the construction of convexity and inflection 
preserving splines 

We observe that results on inflection counts, apart from affecting the analysis of inflection cri- 
teria of splines, have significant effect on the analysis for convexity preservation criteria of 
splines. Convexity preservation criteria requires that under certain conditions projection of a 
curve segment 7(2) on planes with a specified normal N c should be convex. Also the condition 
that inflection count of curve 7(2) is greater than 1, that is, 1(7) > 1 says that there exist a 
vector A^ s , such that, projection of curve on planes with normal vector N s has inflection points 
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greater than 1 and hence is not convex. Thus we see that if N c = N s , then j(t) fails to satisfy 
the convexity criteria. 

Among the results stated below some of them are stated in jH Goodman, 1991] as corollaries 
we state them as theorems because of their relevance to us. 

Theorem 7.1 £5] Goodman, 1991 ] If-y is a curve, as in theorem < \5.2\) . which is not planar, then 

m > i- 

Using the theorem (14.61) we have following results for cubic Bezier curves. 

Theorem 7.2 [5 Goodman, 1991 ] If'-f is a cubic polynomial curve which is not planar, 1(7) = 
2 

We have seen in previous sections that convexity and inflection counts Bezier and B-spline 
curve are related to the convexity and inflection counts of their control polygons. We state few 
results from [5, Goodman, 1991] for inflection count of polygonal arcs. 

Theorem 7.3 [5 Goodman, 1991] If Pq, P n are not coplanar, then I(T) = 1 if and only if 
V\, V n -i lie in order in a plane sector sub-tending an angle < n. 

Theorem 7.4 /[5 Goodman, 1991 ] If n = 3 and Pq, P3 are not coplanar, then I(T) = 1. 

Remark 7.5 The above theorems negates the general perception about the inflection counts 
and convexity of curves and polygonal arcs. 

Remark 7.6 The proofs provided in /[5] Goodman, 1991 ]for the theorems stated above are con- 
structive, that is, plane on which projection of the curves have inflections points are explicitly 
constructed. 

We now illustrate the difficulties in constructing a convexity and inflection preserving curve 
(as indicated in the theorems above) in the following examples. 

Example 7.7 Let^ 2 = (-3,-3,-0.5), Xj_i = (0,0,0), x; = (0,0,5), x m = (2,-4,5.5). 
The normals at Xj_x and Xj are Nj_! = (1.5, —1.5, 0) and Nj = (2., 1., 0) respectively. Since 
iVj_i • Ni = 1.5 > Owe require that the curve between Xj_i and x« satisfy convexity preservation 
criteria. 
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Data polygonal arc D formed by Xj_ 2 , x i-i> x « ««rfx i+1 along with normals Nj_i anJNj as 



?/zzcfc Zmes is shown infigure \13(a)\ The figures \13(b)\ \13(c)\ and \13(d)\ show the curve between 
Xj_x and Xj with data polygonal arc D along the viewpoints Nj_i and Nj and Vi = (—10, 0, 1) 
(that is projection in the plane with normals N^_i, Nj and ViJ. Wfe observe that though the 
curve satisfies the conditions of convexity preservation criteria along the viewpoint Vi, it fails 
to do the same for both the viewpoints N,_i and Nj. TTzws we see it is relatively difficult to 
construct a curve satisfying the convexity preservation criteria using a graphical interface. 



The mathematica code for the generation of figures 13(a) 13(b) 13(c) and 13(c) is as be 



low: 



(*** 



ml={3,l,0.5};m2={2,- 1,0.5}; 

A= {0,0,0}; B={0,0,5}; Am=A+ml; Bm=B-m2; 

P=((l-tr3)*A + (3*(l-tr2*t)*Am + (3*(l-t)*r2)*Bm + (f3)*B 

L0={-3,-3,-0.5}; L1=B-A; L2={2,-4,0.5}; F=B+L2; 

Nl=(l/10)*Cross[(-l)*L0,B] N2=(l/10)*Cross[B,L2] DT=N1.N2 

P0=A+t*L0; Pl=A+t*B; P2=B+t*L2; 

(* L0, A, B,F are data points *) 

Poly=Show[Graphics3D[{Line[{L0,A,B,F}], {Thickness[0.010], Line[{A,Nl}]}, {Thickness[0.010], Line[{B,N2}]}}], 
ViewPoint- >{- 1,0,1}, Boxed- >False]; 

Poly0=Show[Graphics3D[Line[{L0,A,B,F}]], ViewPoint- >{-10,0,l}, Boxed- >False]; 
ln0=ParametricPlot3D[P0,{t,0,l}, Boxed- >False]; 
lnl=ParametricPlot3D[Pl,{t,0,l}, Boxed- >False]; 
ln2=ParametricPlot3D[P2,{t,0,l}, Boxed- >False]; 
Crv0=ParametricPlot3D[P,{t,0,l}, ViewPoint- >N1, Boxed- >False]; 
CrvPoly0=Show[{Crvl,Poly0}, ViewPoint- >N1, Boxed- >False]; 
Crvl=ParametricPlot3D[P, {t,0,l}, ViewPoint- >N2, Boxed- >False]; 
CrvPolyl=Show[{Crvl,PolyO}, ViewPoint- >N2, Boxed- >False]; 
Crv2=ParametricPlot3D[P,{t,0,l}, ViewPoint- >{-10,0,l}, Boxed- >False]; 
CrvPoly2=Show[lnO,lnl ,ln2,Crv2] ; 

Display["d:\Gautam_Viewpoint\view0.png", Poly, "PNG"]; 
Display["d:\Gautam_Viewpoint\viewl.png", CrvPolyO, "PNG"]; 
Display["d:\Gautam_Viewpoint\view2.png", CrvPolyl, "PNG"]; 
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(a) Data polygon arc with normals as thick 
lines 



(b) Viewpoint Ni_i 




Display["d:\Gautam_Viewpoint\view3.png", CrvPoly2, "PNG"]; 



( 



) 



Example 7.8 Here we have Xj_ 2 = (—3, —3, —0.5), Xj_x = (0, 0, 0), Xj = (0, 0, 10), x i+1 = 
(2, —4, 10.5). The normals at Xj_x cmd x, are Nj_i = (30., —30., 0) and Nj = (40., 20., 0) 
respectively. Since A^_i ■ iVj = 600 > we require that the curve between Xj_i and Xj satisfy 
convexity preservation criteria. 

The figures \14(a)\ \14(b)\ \14(c)\ \14( d)(14je)\ and \14(f)\ show the curve between Xj_i and Xj 
wzY/z data polygonal arc along the viewpoints Nj_i, Nj, Vi = (7.458, —1.863, —3.506), V 2 = 
(-17.458,1.863,23.506), V 3 = (-7.458,6.863,33.506) and V 4 = (-7.458,30.863,33.506) 
respectively. Observe that the curve along viewpoints Nj_i, Nj, Vi, V 2 w convex where as 
along viewpoints V 3 and V4 is not convex. This means that if values o/xj_ 2 or Xj + i are altered 
such that Nj_i or Nj changes to V 3 or V 4 wz'ta Nj_i ■ Nj > 0, taen tae cnrve doesn't satisfy 
convexity preservation criteria with respect to changed data polygonal arc. 

In addition to the above observation in example we also note that (figure \14(e)ty along the 
viewpoint V 3 = (—7.458,6.863,33.506) the curve has two inflections as shown in figure [771 
(the curvature of the projection of curve along the viewpoint V 5 changes its sign twice). Thus 
i/Xj_2 or Xj +1 are such that Nj_x or Nj respectively are equal to V 3 , with Nj_i • Nj < then 
the curve would not satisfy inflection preservation criteria. 

The mathematica code for the generation of figures \14(a)\ \14(b)\ \14(c)\ \14(d)\ \14(e)\ \14(f)\ 



is as below: 

m 1 x= 1 ;m 1 y = 1 ;m 1 z=0. 5 ; 
m2x=2 ;m2y=-3 ;m2z=0. 5 ; 
Ax=0;Ay=0;Az=0; 
Dx=0;Dy=0;Dz=10; 

Bx=Ax+mlx; By=Ay+mly;Bz=Az+mlz;(*(l, 1,0.5)*) 
Cx=Dx-m2x;Cy=Dy-m2y; Cz=Dz-m2z;(*(-2,3,9.5)*) 
Px=Ax*(l-tr 3 + Bx*3*(l-t) A 2*t + Cx*3*(l-t)*r 2 + Dx*f 3; 
Py=Ay*(l-tr 3 + By*3*(l-tr 2*t + Cy*3*(l-t)*f 2 + Dy*f 3; 
Pz=Az*(l-tr 3 + Bz*3*(l-tr 2*t + Cz*3*(l-t)*f 2 + Dz*t" 3; 
L0x=-3;L0y=-3;L0z=-0.5;(*xi-2=(-3,-3,-0.5)*) 
L 1 x=Dx- Ax ;L 1 y=Dy- Ay ;L 1 z=Dz- Az ; 
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L2x=2 ;L2y=-4;L2z=0.5 ; 

P0x=Ax+L0x*t;P0y=Ay+L0y*t;P0z=Az+L0z*t; 

Plx=Ax+Dx*t;Ply=Ay+Dy*t;Plz=Az+Dz*t; 

P2x=Dx+L2x*t;P2y=Dy+L2y*t;P2z=Dz+L2z*t;(*xi+l=(2,-4,10.5)*) 
Nl=Cross[{-LOx,-LOy,-LOz}, {0,0,10}] N2=Cross[{0,0,10},{L2x,L2y,10+L2z}] A=N1.N2 
CrvPolyl=ParametricPlot3D[{{Px,Py,Pz}, {P0x,P0y,P0z}, {Plx,Ply,Plz}, {P2x,P2y,P2z}}, {t,0,l}, ViewPoint- 
Boxed— > False]; 

CrvPoly2=ParametricPlot3D[{{Px,Py,Pz}, {P0x,P0y,P0z}, {Plx,Ply,Plz}, {P2x,P2y,P2z}}, {t,0,l}, ViewPoint- 
Boxed— > False]; 

CrvPoly3=ParametricPlot3D[{{Px,Py,Pz}, {P0x,P0y,P0z}, {Plx,Ply,Plz}, {P2x,P2y,P2z}}, {t,0,l}, ViewPoint- 
-1.863, -3.506}, Boxed- >False]; 

CrvPoly4=ParametricPlot3D[{{Px,Py,Pz}, {P0x,P0y,P0z}, {Plx,Ply,Plz}, {P2x,P2y,P2z}}, {t,0,l}, ViewPoint- 
17.458, 1.863, 23.506}, Boxed- >False]; 

CrvPoly5=ParametricPlot3D[{{Px,Py,Pz}, {P0x,P0y,P0z}, {Plx,Ply,Plz}, {P2x,P2y,P2z}}, {t,0,l}, ViewPoint- 
7.458, 6.863, 33.506}, Boxed- >False]; 

CrvPoly6=ParametricPlot3D[{{Px,Py,Pz}, {P0x,P0y,P0z}, {Plx,Ply,Plz}, {P2x,P2y,P2z}}, {t,0,l}, ViewPoint- 

7.458, 30.863, 33.506}, Boxed- >False]; 

Display ["new view 1 .png",CrvPoly 1 ,"PNG"] 

Display["newview2.png",CrvPoly2,"PNG"] 

Display["newview3.png",CrvPoly3,"PNG"] 

Display["newview4.png",CrvPoly4,"PNG"] 

Display["newview5.png",CrvPoly5,"PNG"] 

Display["newview6.png",CrvPoly6,"PNG"] 

The mathematica code for the generation of figure 1731 for the curvature of the projection of the 
curve along the viewpoint V3 is as below: 

ml={l, 1,0.5}; m2={2,-3,0.5}; 

A={0,0,0}; B={0,0,10}; Am=A+ml; Bm=B-m2; 

P=((l-tr 3)*A + (3*(l-tr 2*t)*Am + (3*(l-t)*r 2)*Bm + (C 3)*B 

DRl=D[P,{t,l}] DR2=D[P,{t,2}] CTR=Cross[DRl,DR2] 

CTRPerp=CTR.{-7.458, 6.863, 33.506} 
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(f) Viewpoint V4 = 

(e) Viewpoint V 3 = (-7.458, 6.863, 33.506) 

(-7.458,30.863,33.506) 

Figure 14: Data polygonal arc with curve between Xj_i and along different viewpoints 



Infl5=Plot[CTRPerp,{t,0, 1 }] 

Display["d:\Gautam_Viewpoint\infl5.png", Infl5, "PNG"]; 



( 



) 



Remark 7.9 Note in that in the above examples as we change the viewpoints the curve and 
the data polygonal arc approximately same nature of deviation from convexity. In figures the 
curve as well as data polygonal arc shows approximately same inflections. This due to the fact 
that the curve satisfies torsion preservation criteria described in section along with convexity 
preservation criteria. 



8 Collinearity preservation criteria for interpolating splines 

Definition 8.1 /(9] Karavelas andKaklis, 2000] The collinearity preservation criteria is defined 
by the condition that if\Ni\ = and Lj_i • Li > 0, then 

W{t)xLj\ 



\l'(t)\\Lj 



<e ,ter]i,j = i-l,i, (8.1) 



where e is a user-specified small positive number in (0, 1], and rji a user specified closed subin- 
terval of(ti-i, £i+i) that includes ti as an interior point. 
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Figure 15: Curvature of the projection of the curve along the viewpoint V 3 

We note that |A^| =0 =>- Lj_ x = a G R and this condition with L^i ■ Lj > 
implies that = aLi, a G R + . Equation (18.11 ) states that the (sine of the) angle between 
tangent vector at each point on the spline ^(t) and L is less than e > in the user specified closed 
interval in (tj_i,ij + i). Thus collinearity preservation criteria requires that if two consecutive 
polygon segments are collinear and is having the same direction then the curve segments of 
the corresponding indexes should be approximately collinear and parallel to the corresponding 
polygon segments. 

We now investigate collinearity preservation criteria a bit more closely. One natural question 
to ask is: What if one considers X; to be redundant? Well in that case, if the curve ^(t) is not 
collinear to the line segment {xj_i, x i+ i}, then the curve between Xj_i and x i+ i is required to 
satisfy other shape preservation criteria. For convenience of understanding the situation let us 
suppose Xj_ 2 , Xi_i, Xj, x i+1 , x i+2 are coplanar. Now consider the following cases 

case i A^_i • iVj + i > to be refered as convex neighbourhood data 

case ii iVj_i • N i+ i < to be refrered as inflection neighbourhood data 

For case i we propose that if the curve does not coincide with line segments {x, ; _!, x, }, {x, ; , 
x i+1 } one must ensure that 

1. 7(t) does not interpolate Xj, 

2. uj(tj) ■ Nj > 0, j = i - 1, i + 1, 

3. (V(t i+ i) x L i+1 ) ■ (j'(t i+1 ) x L i+2 ) < and (f^i-i) x U) ■ (f^-i) x L^ x ) < 
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4. j(t) is globally convex between Xj_i and x i+1 , 



5. = aLi, a G R + for t 6 (t* — 77, t« + 5) C (suitable choice of 77 and 5 

provides necessary tilt to the curve 7(t)). 



For a suitable choice of e > condition 18. II along with conditions 1-5 the curve 7(£) have 
following properties (see Figure [Hand [2]): 



satisfies convexity preservation criteria between Xj_i and x i+ i and 

convexity and inflection criteria preservation criteria for data arc segments {x;_i, Xj_ 2 } 
and {x i+ i, Xj +2 } achievable. 



Violation of any of these conditions leads to the violation of the above properties as illus- 
trated by Figure [3] and collinconvwr2. 




X. T 
l+l X 




Figure 16a : Collinearity preservation criteria for convex neighbourhood data makes condi- 
tions for convexity preservation criteria between Xj_i and Xj_ 2 achievable. 




X. T 
l+l X 




Figure 16b : Collinearity preservation criteria for convex neighbourhood data makes condi- 
tions for inflection preservation criteria between Xj_i and x^_ 2 achievable. 
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Figure 16c : Violation of collinearity preservation criteria for convex neighbourhood data 
leads to violation of convexity preservation criteria between Xj_i and x i+1 





Figure 16d : Violation of collinearity preservation criteria for convex neighbourhood data 
leads to violation of convexity preservation criteria between x^_i and x i+ i 

For case ii we propose that if the curve does not coincide with line segments {xj_i, x^}, {x,, 
x i+1 } one must ensure that 

1. j(t) interpolate x,, 

2. uj(tj) ■ Nj > 0, j = i - 1, i + 1, 

3. (Vfti+O x L i+1 ) ■ {i{t i+1 ) x L i+2 ) < and (Y^-i) x U) ■ (YC^-i) x L^) < 

4. uj(t) changes sign only once between Xj_i and x i+ i 

5. u>(t) changes sign at t = tj. 



For a suitable choice of e , condition 18. II along with conditions 1-5 the curve j(t) have 
following properties (see Figure [Hand [2]): 

• satisfies convexity preservation criteria between Xj_i and x i+ i and 
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• convexity and inflection criteria preservation criteria for data arc segments {xj_ l5 x;_ 2 } 
and {x i+1 , x i+2 } achievable. 



Violation of any of these conditions leads to the violation of the above properties as illus- 
trated by [3] and H. 




Figure 17a : Collinearity preservation criteria for inflection neighbourhood data makes con- 
ditions for convexity preservation criteria between x*_i and Xj_ 2 achievable. 




Figure 17b : Collinearity preservation criteria for inflection neighbourhood data makes con- 
ditions for inflection preservation criteria between 3q_i and x^_ 2 achievable. 
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Figure 17c : Violation of collinearity preservation criteria for inflection neighbourhood data 
leads to violation of inflection preservation criteria between Xj_i and x i+i 




Figure 17d : Violation of collinearity preservation criteria for inflection neighbourhood data 
leads to violation of inflection preservation criteria between x^ and x i+1 

Thus we see that we need to modify the definition of collinearity preservation criteria for 
the general data, that is, when Xj_ 2 x i _ 1 , x i5 x i+1 and x i+2 are nonplanar. The modification is 
to be done by adding conditions according to the following cases: 

case i iV7i ■ N i+ i > to be refered as convex neighbourhood data 

case ii Ni_i ■ N i+1 < to be refrered as inflection neighbourhood data 

We state our modified definition as follows: 

Definition 8.2 The collinearity preservation criteria is defined by the condition that if \ Ni \ = 
and • Lj > 0, then 



j(t) x L 



<e ,te rji, j = i - 



(8.2) 



31 



where e is a user-specified small positive number in (0, 1], and rji a user specified closed subin- 
terval o/(ij_i, t i+ i) that includes tj as an interior point and additionally for the case of convex 
neighbourhood data 

1. 7(t) does not interpolate x i; 

2. 7(t) should satisfy convexity preservation criteria between Xj_i and x i+ i, considering 
Xj_ 2 > Xj_x, x i+1 and y. i+ 2 as consecutive data points, 

3. |7 ; (*i+i) x L i+1 \\>y'(t i+1 ) x L i+2 \ < and \i{U-i) x i i ||Y(i i _i) x L;_i| < 

4. j'(t) = aLi, a G R + for t G (t« — 77, + 5) C [it_i,ii+i] (suitable choice of r\ and 5 
provides necessary tilt to the curve j(t) ). 

for the case of inflection neighbourhood data 

1. j(t) interpolate Xj, 

2. 7(t) should satisfy inflection preservation criteria between Xj_i and Xj+i, considering 
Xj_ 2 , Xj_i, x i+ i and x i+2 g« consecutive data points, 

3. |7 ; (*i+i) x L i+1 \\^'(t i+1 ) x L i+2 \ < and \i{ti-i) x LiHY^i) x £i_ a | < 

Considering Xj_2, Xj_i, Xj+i and Xj + 2 a* consecutive data points the curve j(t) between Xj_i 
and Xj+i a/50 satisfy torsion preservation criteria ( to be stated in section [9]) or coplanarity 
preservation criteria (to be stated in section 1701) according to the condition [Lj_i Lj Lj+i] 7^ 
or [Lj_i Lj = respectively. 

We observe that the last two conditions guides the spatial behavior of the curve 7(f) be- 
tween Xj_x and x i+1 , with respect to the planes IIj_i (containing Xj_ 2 > Xi-i> and Xj) and 
(containing Xj, x i+ i and x i+2 )- We also observe that in case x, ; _ 2 , Xj-i, Xj, Xj +1 and x i+2 are 
nonplanar, then conditions of the definition 18.21 makes conditions of coplanarity preservation 
criteria (to be stated in section [10]) for curve 7(2) between x i+2 and x i+3 and between Xj_ 2 and 
Xj_i, achievable. Thus we see that for a person, who tends to ignore Xj as a data point and 
considers Xj_i and x i+1 as adjacent data points, definition 18.21 makes all the shape preservation 
criteria by j(t) between Xj_i and x, + i achievable, without conflict with the shape preservation 
criteria for curve segment between Xj_ 2 and Xj_i and between x i+ i and x i+2 . 
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9 Torsion preservation criteria for interpolating splines 



Definition 9.1 Discrete torsion for the polygonal arc X0X1 • • • x n is defined as 

Aj = L, L i+ i], z = 3, . . . , n - 1 (9.1) 
where Lj = Xj — X£_i, 2 = 1, n 

Definition 9.2 / f7?l Costantini and Manni, 2003] Torsion preservation criteria consists of fol- 
lowing conditions 

1. Ti(t) Aj > in a chosen closed subinterval of(ti~i,ti), whenever Aj 7^ 0. 

2. Tj(^_x)Aj > 0, j = i — 1, i, whenever Aj_xAj > 

where Ti{t) = l W^iwF ' ¥im x J ' m + °" 

First condition of torsion preservation criteria states that i th curve segment should appear 
to twist away from its osculating plane in the same way as L i+ i moves away from the plane of 

{Xj_ 2 , Xj_i, Xj}. 

Most of the author in their papers do not consider second condition in their definition for 
torsion preservation criteria. Following definition is followed by them 

Definition 9.3 / fiJ] Kong and Ong, 2002] Torsion preservation criteria is defined by the con- 

W.(t) <y" (t) 1"' (t)\ 

dition that if A, ^ then r i (t)A i > 0, t e [U-iM' where n{t) = \, ' „L ui > tf 

iS) x i'm * 0. 

(We discuss the situation in which Tj(t) Aj = in the theorem [T2J]) 

We now discuss the second condition of definition of !9.2[ In figure [T8l set of points with bigger 
circle, {xj_ 2 > Xj_i, Xj, x i+ i}, correspond to Aj and set of points with smaller circle, {xj_ 3 , 
Xj_ 2 , Xj_i, Xj}, correspond to Aj_i. We know that the sign of Aj and Aj_i depends on the 
cosine of the angle that L i+ \ and Lj_ 2 , respectively, makes with the normal of the plane Ilj_i 
containing Lj_i and Lj. The condition AjAj_ x > states that Lj_ 2 moves into the plane Ilj_i 
in the same way as L i+1 moves out of the plane, that is, the line segments {xj_ 3 , Xj_ 2 } and {xj, 
x i+1 } lies on the opposite sides of the plane ITj_i. The condition Tj(tj_i) Aj > 0, states that the 
curve at t = moves out of its osculating plane in the same way as the vector L i+ i moves out 
of the plane ilj_x. Similarly, the condition Tj(tj_i)Aj_i > 0, states that the curve at t = tj_i 
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Figure 18: Data polygon arc involved in second condition of definition 19 .21 



moves out of its osculating plane in the same way as the vector Lj_ 2 moves into the plane 
or Li moves out of the plane Hi-2- 

We first observe that sign(rj(t i _ 1 )Aj_i) = sign(r i (tj_ 1 )A i ), whenever Aj_iAj > and 
sign(T i (4_i)A i _i) = -sign(r i (t i _ 1 )A i ), whenever A^A* < (since sign(r i (t i _ 1 )A i _ 1 ) = 
sign^^^A^A,) 2 ) = sign(r i (ti_i)A i (A i _ 1 A i )) ). 

Given Aj_i Aj > 0, the condition Ti(tj_i) Aj_i > implies that Tj(tj_i) Aj > and vice-versa. 
The conditions r^i^^Aj > and Aj_]Aj < implies that rj(t i _ 1 )A i < 0. For the case 
Aj.xAj < 0, that is, the line segments {x.;_ 3 , x.j_ 2 } and {xj, x i+1 } lying on the same side of 
the plane Hi-i, the curve satisfying the condition 7i(tj_i) Aj > 0, have the property that 

• (since condition rj(£j_i)Aj > 0) the curve at t = moves out of its osculating plane in the 
same way as the vector L i+1 moves out of the plane ITj_i. 

• (since rj(t i _ 1 )A i _ 1 < 0) the curve at t = moves out of its osculating plane opposite to 
the way as the vector Lj_ 2 moves into the plane or Li moves out of the plane ilj_2. 

Since in the definition [93j we have Tj(tj) Aj > 0, above analysis holds true when the curve 
is traversed in the reverse direction by concentrating the view on %i instead of Xj_i (with set of 
points involved being {xj_ 2 , Xi-i, x, ; , x m , x i+2 }) 

Thus from the above analysis, we see that the definition [9J] has better influence on the shape 
of the curve than that of the definition 19 .21 However, the definition 19 . 21 helped us to bring a very 
important property (discussed above) of curve satisfying condition of definition 19.31 
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10 Coplanarity preservation criteria for interpolating spline 



Definition 10.1 [9: Karavelas and Kaklis, 2000] Coplanarity preservation criteria is defined 
by the condition that if Aj = and |A^_i||A^| ^ 0, then 



where e\ is a user specified small positive number in (0, 1], and U is user-specified closed 
interval such that C I i C (ij_2, t i+ \). 

Coplanarity preservation criteria states that if data points Xj_ 2 , x i-i> x i an d x i+1 are coplanar 
to a plane II, then the interpolating curve between Xj_i and x, and in the vicinity of Xj_i and Xj 
has its binormal close to Nj, that is, its osculating plane should remain close to a plane parallel 
to U r 

We observe that in addtion to the condition (110.11) . the curve segment between Xj_ x and Xj 
should be constrained such that its oscillations about the plane EL, is minimum. In fact, 

• if \u(ti-i) x Ni-i\ = and \u(ti) x iVj| = 0, then the curve segment can be constrained to 
be coplanar with the plane Ilj. 

• if x Ni-i\ ^ and \u(U) x Ni\ = (or |o;(ti_i) x jV<_i| = and \uj(U) x ^ 0) 
then the curve segment can be constrained such that it oscillation about the plane Ilj only once. 

and if |o;(tj_i) x N^ x \ ^ and \uj(U) x ^ 

• and addtionally iVj_i and N{ lie on the same side of the plane Ilj, then the curve segment can 
be constrained such that it does not oscillate about a (fixed) plane parallel to the plane Ilj and 

• and addtionally when iVj_i and N{ lie on the opposite side of the plane Ilj, the curve segment 
can be constrained such that it oscillates about a (fixed) plane parallel to the plane Ilj only once. 

11 Different shape preservation criteria on a curve segment 

We now observe that the data points {xj_ 2 , Xj_i, Xj, Xj + i} must satisfy one of the two conditions 
cl iVj.x • Ni < 0, qualifying condition for inflection preservation criteria, 
c2 iVj_! • Ni > 0, qualifying condition for convexity preservation criteria 
with one of the two conditions 

tl Aj = 0, qualifying condition for torsion preservation criteria, 



Uj(t) X Nj 

Ht)m\ 



< e lt \uj{t)\ ^ 0, t G U j = i - l,i 



(10.1) 
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t2 Aj 7^ 0, qualifying condition for coplanarity preservation criteria. 



We observe that there is no conflict between the conditions that the curve needs to satisfy for 
one among convexity preservation criteria and inflection preservation criteria simultaneously 
with one among torsion preservation criteria and coplanarity preservation criteria. 

We also observe that there is no conflict between condition that curve need to satisfy for 
collinearity preservation criteria simultaneously with torsion preservation criteria or coplanarity 
preservation criteria. 

12 Avoiding conflict between shape preservation behaviour 
of adjacent curve segments 

We observed in previous section that every curve segment need to satisfy 1) condition for either 
convexity preservation criteria or inflection preservation criteria 2) condition for either torsion 
preservation criteria or coplanarity preservation criteria. In this section we investigate the com- 
patibility between the shape preservation behaviour of adjacent curve segments. In Figure 19a- 
19d we observe that if the curve j(t) is required to be C 1 smooth then there is possibility that 
convexity as well as inflection preservation of a curve segment may lead to the violation of 
convexity and inflection preservation of adjacent curve segment. 



Figure 19a : Convexity preserving curve segment between Xj_i and Xj making violation of 
inflection preservation criteria between Xj_i and Xj_ 2 imminent. 




i-3 
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X 

i+1 

Figure 19b : Convexity preserving curve segment between Xj_i and Xj making violation of 
inflection preservation criteria between Xj_x and x,_ 2 imminent. 

X ✓ 




Figure 19c : Inflection preserving curve segment between Xj_i and x 4 making violation of 
inflection preservation criteria between Xj and x i+1 imminent. 




Figure 19d : Inflection preserving curve segment between Xj_x and x* making violation of 
convexity preservation criteria between Xj and x i+1 imminent. 

We observe from Figure 19a-19d that if Xj_ 2 > x i-i> x i> x ;+i are coplanar, then the conflict 
of convexity preservation criteria or inflection preservation of a curve segment with that of ad- 
jacent curve segment (of C 1 smooth spline curve 7(i)) is resolved if and only if (7'(£j_i) x Li) ■ 



37 



(Y(^-i) x L*_i) < 0. 



For the general case, we observe from the definition of convexity and inflection preservation 
criteria that compatibility between convexity and inflection preservation behaviour of adjacent 
curve segments of C 1 smooth spline curve ^(t) can be guaranteed if and only if 

(l' N ± (U-i) x L { ) ■ ( 7 ^ ± (ti-i) x < (12.1) 

i— 1 i — 1 

where 7^(0 is the orthogonal projection the curve 7(t) on a plane having normal vector as 
iVj. We note that the condition (112.11) does not interfere with the conditions of convexity and 
inflection preservation criteria. 




Figure 20a : Convexity preserving curve segment between x$_i and Xj satsifying (112.11) facil- 
itates inflection preservation between x^ and Xj_ 2 - 




Figure 20b : Convexity preserving curve segment between Xj_i and x, ; satsifying (112.11) facil- 
itates inflection preservation between Xj_i and X;_ 2 . 
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Figure 20c : Inflection preserving curve segment between Xj_i and x, satsifying (112.11) facili- 
tates inflection preservation between x, and x i+1 . 




Figure 20d : Inflection preserving curve segment between Xj_i and x, satsifying (|12.1I) facil- 
itates convexity preservation between x» and Xj + i. 

We also observe that the condition (112.11) can be conveniently imposed on a curve along 
with the conditions of torsion and coplanarity preservation criteria. Therefore, compatibility 
of torsion and coplanarity preservation of a curve segment with the convexity and inflection 
preservation of adjacent curve segment is feasible. From the above analysis and conditions of 
collinearity preservation criteria, we see that compatibility of convexity, inflection, torsion and 
coplanarity preservation of a curve segment with the collinearity preservation of adjacent curve 
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segment is also feasible. 

We now analyze the compatibility between the adjacent curve segments satisfying copla- 
narity and torsion preservation criteria. Let us consider figure [T8] The points Xj_ 3 , x.j_ 2 , Xj-i, 
Xj and Xj+i may be such that 1) AjAj_i < 0, 2) AjA,;_i > and 3) A, ; ^ 0, Aj_i = 0. 

Theorem 12.1 If A^Aj_i < 0, then the curve segments 7 j_i(t) (between Xj_ 2 and ^i-i) and 
7i(i) (between x$_i and xj, of spline curve j(t), satisfies conditions for torsion preserva- 
tion criteria, if and only if either j(t) is torsion discontinuous or r(^_i) = (and therefore 
niU-^Ai = 0, ri-i^-OAi-! = 0). 

Proof: As discussed in section [9] we have sign(7i(tj_i)Aj) = sign(7i(£j_i) Aj_iAjAj_i). 
Therefore, from the definition |9.3| we see that torsion preservation by 7j(t) requires Tj(tj_i) A^ < 
and torsion preservation by 7i_i(t) requires Tj_i(tj_i) Aj_i > 0. Hence the theorem. I 
From the above proof it is evident that if A, Aj_i > 0, then torsion preservation by 7i _i and 
ji is compatible. 

If Aj 7^ and Aj_i = (requiring coplanarity preservation by 7i_i(i)), then from the 
definition 1 10. II we see that 7 j(t) need to satisfy the condition of coplanarity condition llO.li that 
is, binormal of the curve should be close to AT i _ 1 in addition to satisfying torsion preservation 
criteria, fort G [U-^ti] H k-u where [U- 2 ,U-i] Q h-i C (U_ 3 ,ti). 

13 Shape preservation by cubic interpolating splines 
13.1 Cubic Bezier segments 

Let the control polygon for Bezier representation of i th cubic curve segment ^(t) of cubic 
spline be {Pj.o, Pj,i, Pj,2> Pi.3}> with P i = Xj and P i 3 = x i+1 . For the i th cubic curve 
segment ji(t), of the cubic spline 7 (t) we have 

7i (t) = P ij0 B*(u(t)) + P <fl S3(«(t)) + P h2 B 3 2 (u(t)) + P i>3 Bl(u(t)) (13.1) 

where B'^{t) t G [0, 1] is 2 th Bernstein's polynomial of order n, = Conditions for 

shape preservation criteria consists of first, second and third order derivatives of Bezier curves. 
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Therefore we get their expressions in terms of the end point and slopes at end points. 

i-Xt) = ^((P;,! - P ifi )B 2 (u(t)) + (P <)2 - P itl )B*(u{t)) + (P il3 - P l , 2 )B 2 2 (u(t))) 

(13.2) 

7 f(t) = ^((p <t2 -2P i) i + P <i o)(l-u(t)) + (Pi,3-2P iia + P ii i)u(*)) (13.3) 
7 f(t) = ^(P^-3P <i2 + 3P iil -P i , ) (13.4) 

itt , 3(Pj i — Pj o) 3(Pj 3 — Pj 2 ) , r t> T-u f 

We have mj_i = ; , rrij = — ■ — and Lj = Pj 3 — Pj . Therefore we can 

hi hi 

rewrite the expression for curve and its derivatives as follows 

7i(t) = x^B^uit)) + (xi_ a + ^mi_i)B?(«(t)) + 

(Xj - y mi )5 2 3 ( M (t)) + x.Sf (>(t)), (13.5) 

7 '(t) = m^^lXt)) + (-A - m^x - mi)Bl{u(t)) + m^ 2 ^)), (13.6) 

l"{t) = ^((-L < -2m < _ 1 -m i ))(l-tz(t)) + (--L i + m i _ 1 + 2m j )u(*)) 



(13.7) 
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^"(*) = T3(/K(m<-i + m<)-2Li). (13.8) 



13.2 Convexity preservation criteria for cubic interpolating splines 

Recall from theorem ?? that cubic spline 7 (i) interpolating data points x, ; , < i < n, satisfies 
convexity criteria if and only if the control polygons of the projection of 7 (t), t G tj] on 
planes with normal vectors Nj, j = i — 1, % are globally convex, whenever iVj_i • iVj > 0. 
We now find simplification of the condition of global convexity control polygon of P/v±( 7 (t)), 
£ G Using this simplified condition we find the modified convexity preservation criteria 

for cubic splines in theorem fl 3 .4[ 

Lemma 13.1 [3, Goldman, 1990] Let the points Pq, Pi P 2 , P3 in R 3 lie on plane with normal 
vector n. Then a line, through the points Po and Pi, intersects with a line, through the points 
P 2 and P3, at the point 

P= P +(Pi-P )s= P 3 + (P 2 -P 3 )t 
= P 1 + (P -Pi)s= P 2 + (P 3 -P 2 )I 
(P 3 - P ) x (P 2 - P 3 ) • n (P x - P ) x (P 3 - P ) ■ n 

~9 v 



(Pi - P ) x (P 2 - P 3 ) • n ((Pi - P ) x (P 2 - P 3 )) • n 
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(P 2 - Pi) x (P 3 - P 2 ) ■ n (P - PO x (P 2 - PO • n 



(P - Pi) x (P 3 - P 2 ) • n (P - Pi) x (P 3 - P 2 ) • n 

Proof: To find the point of intersection the given two lines we need to solve the equations 

P = P + (Pi-P )s (13.1) 
P = P 3 + (P 2 -P 3 )t (13.2) 

for s and t with the condition for the coplanarity of the four points 

(P 3 -Po)-((Pi-Po)x(P 2 -P 3 )) = 

By subtracting equation (113.11) from equation (113.21) we get 

(P 3 -P ) + (P 2 -P 3 )t-(Pi-Po)s = (13.3) 

Taking the cross product on both sides of equation (113.31 ) by (P 2 — P 3 ) we get 

(P 3 - P ) x (P 2 - P 3 ) = ((Pi-Po) x(P 2 -P 3 ))s (13.4) 

Now taking the scalar product on both sides of equation (113.41) with n we get 

((P 3 -P )x(P 2 -P 3 ))-n 

((Pi-P )x(P 2 -P 3 ))-n 1 ^ J 

Similarly we get the value for t. Now by interchanging P with Pi and P 3 with P 2 we get the 
values for s and t. Hence proved. I 

Lemma 13.2 A planar polygonal arc {Po, Pi, P 2 , P 3 } lying on a plane with normal vector N 
is globally convex according to orientation induced by N if and only if either 

1. (Pi - P ) x (P 2 — P 3 ) • iV > with 

(a) (Pi - P ) x (P 2 — Pi) ■ N < and (P 2 - P x ) x (P 3 - P 2 ) • N < or 

(b) (P - Pi) x (P 3 - P ) ■ N < and (P 3 - P ) x (P 2 - P 3 ) ■ N < 

or 

2. (Pi - P ) x (P 2 - P 3 ) ■ iV < with 

(a) (Pi - P ) x (P 2 - P x ) • N > and (P 2 - P x ) x (P 3 - P 2 ) • iV > or 

(b) (P - Pi) x (P 3 - P ) • N > and (P 3 - P ) x (P 2 - P 3 ) • N > 
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holds. 



Proof: From the definition 13 .41 we know that planar polygonal arc {Po, Pi, P2, P3} is globally 
convex according to the orientation induced by normal vector N if and only if 

condition i polygonal arc starting from P always turn towards the right side and 

condition ii it always lies entirely to its right side of any of its edges. 

For the given polygonal arc, condition i holds if and only if 

((Pi - P ) x (P 2 - PO) ■ N((P 2 - Pi) x (P 3 - P 2 )) ■ N > (13.6) 

Now for the given polygonal arc, condition ii holds if and only if 

1. line through P and Pi does not intersect the line segment between P 2 and P 3 and 

2. line through P 2 and P 3 does not intersect the line segment between P 1 and P . 

The above two condition holds if and only if the point of intersection P, between the line l\ 
through P and Pi and line Z 2 through P 2 and P3, does not lie in the segment between P and 
Pi or the segment between P 2 and P 3 . From lemma [TXTl we know that point of intersection of 
lines li and l 2 is given by 

P= P + (Pi-P )s= P 3 + (P 2 -P 3 )t 
= P 1 + (P -P 1 )s= P 2 + (P 3 -P 2 )t 

= ((P3-P0) x (P 2 -P 3 ))-N t = ((Pi - Pp) x (P 3 - P )) ■ N 
((P1-P0) x (P 2 -P 3 ))-N' ((P1-P0) x (P 2 -P 3 ))-N 

_ = ((P2-P1) x (P 3 -P 2 ))-N = ((P0-P1) x (Pa-P^-N 
3 ((P0-P1) x (P 3 -P 2 ))-N' ~" ((P0-P1) x (P 3 -P 2 ))-N 

Now with condition i ensured, condition ii is satisfied if and only if either s < with t < 
or s < with t < holds. Hence proved. I 

Lemma 13.3 Let Pi = P Jv j_(P i ). Then 

(P -P 6 ) x (P c -P d )-iV = (P a - P 6 ) x (P c - P d ) • N (13.7) 
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(a) case 1 (a) s < 0, t < 



(b) case 1 (b) s < 0, t < 







\ P 



\ ' 
<P 
\ 



(c) case 2 (a) s < 0, t < (d) case 2 (b) s < 0, t < 

Figure 21: Examples of different planar convex control polygons 



Proof: We know that projection of Pj on a plane 



(x, y,z) ■ N + d 



Pi is given by 



Pi + p -; N J d N . 



with normal vector N, 



(13.8) 



One can get the proof using the idea in the proof of lemma I3T31 I 
Using the definition 13.161 and following theorem from 0TJ Liu, 2001] we get 

Theorem 13.4 A cubic spline curve j(t) satisfies the convexity preservation criteria if and only 
if either 

1. rrij x x nij ■ Nj < 0, with 

hi hi 

(a) rrij_i x L, • N« < — nij_i x rrij ■ Nj and Li x rrij ■ N < — nij_i x mj ■ Nj or 

3 3 

(b) mj_i x Li ■ Nj > and Li x rrij • N > 



or 
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2. nij_i X nij • Nj > 0, with 

hi hi 

(a) mi i x Li ■ Nj > — m^i x m, ■ Nj and Li x nij ■ Nj > — m 4 _! x m, • Nj or 

3 3 

(b) m^x x Lj ■ Nj < and Lj x rrij ■ Nj < 
j — i — whenever N^-x ■ Ni > 0. 

Proof: We note that 

(Pi - P ) x (P 2 — Pi) • N — (Pi - P ) x (P 3 - P ) • N + (Pi - P ) x (P 2 - P 3 ) • N 
(P 2 - Pi) x (P 3 - P 2 ) • JV = (P 3 - P ) x (P 3 - P 2 ) • N - (Pi - P ) x (P 3 - P 2 ) • JV 

since P 2 - Pi = (P 2 - P 3 ) - (Pi - Po) + (Ps - Po) 

Since control polygon of (orthogonal) projection of Bezier curve on a plane is same as the 
orthogonal projection of control polygon of the curve on the plane, using theorem ??, lemma 
113.21 (by replacing Pj by Pjj) and lemma [T331 we get the theorem. I 

13.3 Inflection preservation criteria for cubic interpolating splines 

The definition 16 .21 of inflection preservation criteria involves the curvature term Ui(t). So, in or- 
der to get a simplified characterization for inflection criteria for cubic case we first get simplified 
expression for uJi{t) as follows. 

Lemma 13.5 Let c(t) = P (l — t) 2 + Pi (2 t (1 — t)) + P 2 t 2 be a quadratic Bezier curve. Then 

c(t) x c'{t) = 2(P x P x )(l - 1) 2 + (P x P 2 )(2 t (1 - t)) + 2(Pi x P 2 )t 2 . (13.1) 

Proof: In the statement of the lemma we observe that though the c(t) and c'(£) are Bezier curves 
of degree 2 and 1 respectively, their cross product is a Bezier curve of degree 2 instead of 3. In 
order to understand this we express the curve c(t) in power basis form as c(t) = po + pit + P2^ 2 , 
where po = P , Pi = 2(Pi-P ),p 2 = P -2Pi + P 2 . Now the z-coordinate of (c(t) x c'(t)) 
is 



P0,x + Pl,xt + P2,xt 2 P0,y + Pl,yt + p 2 ,yt 2 



Pl,x + 2P2,xt 



Pl,y + 2P2,yt 





Pl,x 


Pi,y 


t 2 + 2 


P0,x P0,y 




P2,x 


P2,y 




P2,x P2,y 



P0,x P0,y 
Pl,x Pl,y 
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The coefficient of t 3 is due to special relation between the curve and its derivative. (We study 
this phenomenon, in detail, in chapter ??.) Thus we see that c(i) x c'(t) = (p x x P2)t 2 + (po x 
P2)t + (Po x Pi)- On substituting the values of pi we get the relation [TXTl I 
Now by substituting the expression for ^[{t) from section [T3TT1 we get the expression for 
curvature of cubic curve as 

<t) = 7K*)X7^) 

= gb,i(l - u(t)) 2 + gi (i «(t)(l - «(*)) + g2,Mt)) 2 (13.2) 



where 



6 2 

g ,i = J^( m i-l X L i) - /H 111 *-! X m *) ( 13 - 3 ^ 

2 

gii = v-Cm^i x mi) (13.4) 



6 2 

g2,i = ^(^i x mi) - ^-( m i-i x m i) (13.5) 



Theorem 13.6 If the i th curve segment of interpolating spline 7$(t) is a cubic curve, U)i(tj) ■ 
Nj > and cjj(t) ■ Nj > changes sign only once for t G [V-i, J = i — 1, i then for all 
N = A./Vj_i + /iA^j, where \fi < 0, cjj(t) ■ iV /las precisely one sign change for t G £»]. 

Proof: We prove the theorem for the case A > 0. We first note that AcUj(tj_i) • Aj_i > 
/MjJi(ti-i) ■ Ni > 0. Since ^(t) ■ Nj > changes sign only once for £ G U], j = % — 1, % 
we also have Ao^tj) ■ A^-i < [iWiiti) ■ N < 0. Thus Aa»j(ij_i) • A^_i + /ia;j(^_i) • iVj = 
• N > and A^(^) • iVj_i + /i^(^) • = u^) • AT < 0. 

From formula (113.21) we see that uJi(t) is a quadratic curve and hence u>i(t) ■ N is quadratic 
polynomial for j — i — 1, i. That is, cjj(t) ■ A" can change sign only twice. And if U{(t) ■ AV-i 
changes sign twice for £ G [tj-i, then we must have a>j(tj_i) ■ A^ > and ^(tj) ■ AT > 0. 
Therefore Ui(t) ■ Nj changes sign only once for t G [V-i, U). 

The proof for the case A < is similar. I 

Remark 13.7 In [\7< Goodman and Ong, CAGD 15, 1-17, 1997], theorem Yl3.6\ is proved for 
rational cubic Bezier curve for a special case in which the tangent vector uij lies on the plane 
with normal vector as Nj, for j — % — 1 , i. But this is a heavy restriction for modeling curves 
and surfaces. 

Theorem 13.8 Ij "the i th curve segment oj "interpolating spline 7«(t) is a cubic curve and Ui(ti-i)- 
Aj_i > 0, andojiifi) ■ Aj_i < then Ui(t) ■ AV-i changes sign only once for t G [U-i, U]. 
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Proof: From formula (113.21 ) we see that u>i(t) is a quadratic curve and hence u>i(t) ■ Nj is 
quadratic polynomial for j — i — 1, i. That is, Ui(t) ■ iV$_i (and u>i(t) ■ Ni) can change sign 
only twice. And if Ui(t) ■ iVj_i changes sign twice for t E [ti-i,U], then ufoi-i) ■ iVj_i and 
uJi(U) ■ A^j_i will have same sign. But we have Uifa-i) ■ iVj_i > and Ui(U) ■ A^_i < 0. 
Therefore uJi(t) ■ Nj changes sign only once for t E i*]. I 

Theorem 13.9 Ifthei th curve segment of interpolating spline 7i(t) is a cubic curve and uJiit^x)- 
Ni < 0, and u>i(ti) ■ iVj_i > then u>i(t) ■ Ni changes sign only once for t E [<i_i, t{]. 

Proof: The proof is similar to that of theorem fl 3 .81 I 

Remark 13.10 The ease of using cubic curve, justified by theorems \13.6\ 17X81 and 17379] is one 

of the reasons for using cubic curve instead of higher order curves in splines. 

Theorem 13.11 ji(t) satisfies inflection criteria if and only if 

1. g ,i ■ iVi_i > 0, g ,j ■ Ni < and 

2. g 2 ,i • < g 2)i • Ni > 
whenever iVj_i ■ Ni < 0. 

Proof: Proof follows from the conditions in the definition 16 .21 for inflection criteria of splines, 
equation (THO) and theorem \TJM fITgl and \TJM I 

Remark 13.12 It is known that quadratic curves cannot be used in splines interpolating non- 
planar set of data points as it does not exhibit torsion. But now we can easily see that the 
quadratic curves cannot even be used in splines interpolating planar set of data points because 
it cannot satisfy the conditions for inflection preserving criteria. 

13.4 Torsion preservation criteria for cubic interpolating splines 

We first state our result for torsion of a cubic Bezier curve. Using the identities in section [TX31 
we get 

Theorem 13.13 Let 7$) = H{t)li{t)i"{t)\, t E fc-i,^] (numerator of r^t)). Thenr^t) = 

T7( m i-1 x m i ■ L i) 
h i 
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Proof: We can write numerator of Tj(t) as 



m = mw(thr(t)\ 

= ui{t)-i!'{t) 

= h ,i(l - u(t)) 2 + h liiU (t)(l - u(t)) + h 2 ,i(u(t)f 



where 



6 2 6 

h,i = (^(mi-i x Li) - ^-( m *-i x m i)) ' 73(^(^-1 + mi) - 2Li) 

6 12 
= ^ i (6(m i _ 1 x L^ ■ uii + 4(mi_i x m*) • L-) = (m i _ 1 x L t ■ mi) 

2 6 24 

h hi = j^( m i-i x m • -^(h^m^x + raj) - 2Lj) = -^(m^x x Lj • mi) 

6 2 6 

^2,i = (^(i, x mi) - -(m^i x mi)) • ^(/i^m,.! + m;) - 2L ; ) 

6 12 
= {6(Li x m^ • mi_i + 4(m i _ 1 x mi) • Li) = j^( m i-i x L i ' m i) 

u(t) = t G \ti-i, U]. Therefore 

ti H—l 

Ti{t) = ^((m*-i x L t - mi )((l~u(t)) 2 + 2u(t)(l-u(t)) + (u(t)) 2 ) 

= ^(irii-i x Li ■ mi) 

Hence proved. I 
From definition [93] of torsion preservation criteria and theorem fl 3 . 131 we get the following 

Theorem 13.14 A cubic spline curve j(t) satisfies torsion preservation criteria, fort G U], 
if and only if [m^ Li mjJAi > 0, whenever Aj 7^ 0. 



13.5 Collinearity preservation criteria for cubic interpolating spline 

With the interpretation given in section [8] we get a sufficient condition for collinearity preser- 
vation criteria for splines in terms of Bezier control points of curve segments as below. Let 
the control points of curve segment j(t), t G [ti, t i+ i] be Pij, j = 0, 1, . . . , rrii and the con- 
trol polygon be represented as Ti. The control points of the derivative of curve segment are 

3 

Pi = — (Pi j + i — Pi A j = 0, . . . , mi — 1 and the control polygon be represented as V[. 

hi 
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Now we make a observation about Bezier curves based on following interpretation about vector 
product between two vectors in R 3 . 

|A x B| 



sin(0), A, B e R , is the angle between A and B (13.1) 



|A||B, 

For any vector L E R 3 equation [13.11 implies the following 

Lemma 13.15 IfOi — Zp*L < 90°, i — 1, 2, p 1; p 2 e -R 3 , then 

[ |p x L| . . r A f Ipi x L\ |p 2 x LI 

SU P \ i I. ri : P = Ph 1 _ + P 2 *' * G i ' !J r = SU P 



IpII-^I " J I IpiII-^I IpsIIl 

Following lemma follows from lemma [T3.15l 



X .£/ 

SU P ^ 1 — uTT '■ P zs a point inside the planar triangle formed by the points 



Lemma 13.16 lfQ { = Zp^L < 90°, t = 1,2, 3, r/zen 
pxl 

IpP 

si f I Pi X L\ 

pi, p 2 and p 3 e i? > = sup < — — — .'i = l,2,3 
We observe that 

Theorem 13.17 For a Bezier curve c(t) with control points p^ 9i = ZpiL < 90°, i 
0,1, ... ,m. 



f left) X LI r A 

p ^ 

SU P "1 — rrrr • P belongs to set of vertices of convex hull formed by Pii = 0, 1, . . . , m 



Pi L 



(13.2) 



Proof: Proof follows from equation (113.11) convex hull property of Bezier curves (which states 
that Bezier curves lie inside the convex hull, that is, smallest convex polyhedra, with triangular 
sides, formed by its control points) and lemma [T3.16l I 

Theorem 13.18 j(t) satisfies collinearity criteria if 

^ : k = 0,1, 2 I < e, j =i - (13.3) 



sup 



\Pi,k\\Lj\ 



whenever \N,j\ = and L;_i • L; > 0, < 9 k — Z.P' i<k Lj < 90°, k — 0, 1, 2, j — i - 1, i (a 
reasonable assumption to make ). 
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Proof: Since for a quadratic Bezier curve convex hull of its control points is the triangle formed 
by the control points, therefore 

{Wit) x LA . A 



sup <{ ' *f : A: = 0, 1, 2 J. , j = % - 1, i, (13.4) 

Hence proved. I 
In case in the definition 1 10. II the condition 1 10. II is replaced by Ui(t) = 0, t E r\i as stated in 
papers Q, Goodman and Ong, CAGD 15, 1997] etc., the collinearity condition for cubic spline 
would have required rrij = aLi, a > in place of 113.31 



13.6 Coplanarity preservation criteria for interpolating cubic splines 

From the definition of coplanarity criteria 110.11 and the analysis in the previous section we 
have following theorem stating the sufficient condition for cubic spline to satisfy co-planarity 
preservation criterion. 

Theorem 13.19 ^(t) satisfies the co-planarity preservation criteria if 

SU P \ I mam :k = 0,l,2\<e,j=i-l,i, (13.1) 

whenever A* = and |JV i _i||JV i | ^ 6 k = Zg itk Nj < 90°, k = 0, 1, 2, j = i - 1 or i (a 
reasonable assumption to make). 

Proof: Proof is same as the proof of theorem [1338] I 
In case in the definition 1 10. II the condition 110.11 is replaced by Ti(t) = 0, t E [U-i,U] as 
stated in papers [7, Goodman and Ong, CAGD 15, 1997] etc., the coplanarity condition for 
cubic spline would have required nij_i x Lj • rrij = 0, that is, nij = aL { + /3mj_x, a, (3 E R in 

place of mm 

Alternatively (actually more precisely) the condition 113.11 can be replaced by the pair of 
conditions mj_i = «iA + /3iLj_i, m 8 = a 2 Li + /5 2 L i+1 , with ol\, a 2 , Pi, P2 > 0. 



14 Shape preserving properties of cubic Catmull-Rom splines 

Let denote the vector x i+ i — Xj_i as tj, i — 2, n — 1, the plane containing the data points 

Xj_i, Xj, Xj + i as Ilj, i — 2, n — 1 and for a curve j(t) = [x(t), y(t), z(t)], t E [0, 1] in R 3 
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let u(t) = j'(t) x 7"(t). We now prove that cubic Catmull-Rom splines, with magnitude of 
tangent vectors calculated according to our algorithm, preserve convexity, inflection, torsion, 
collinearity and coplanarity behavior of the data polygon, as follows. 

The tangent vector of the Catmull-Rom spline at a data point Xj, i = 2, n — 1 is parallel 
to tj so that the tangent vector is coplanar with the plane Ilj. Due to this cubic segment of 
the spline, between the data points Xj and x i+ i, at Xj lies on Ilj and at x i+ i lies on the plane 
Thus from the interpretations torsion preservation and coplanarity preservation criteria 
in sections [9] and [10] we observe that the Catmull-Rom splines (having tangent vectors with our 
magnitudes) preserve the torsion and coplanarity behavior of the data polygon. 

14.1 Torsion and coplanarity preservation 

Torsion preservation by cubic Catmull-Rom spline is also assured by the theorem 113.141 as 
follows. According to the theorem 113.141 the spline curve 7(i) must satisfy the condition 
[m;_i Li rrij]Aj > 0, where Aj = [L»_i Li L i+ i\. But for Catmull-Rom spline we have 
rrij = x i+ i - Xi_i = Li + L i+ i, so that [m;_i Li m*] = [L^i L { L i+1 ] = A» and therefore the 
torsion preservation condition [mj_x L, ; m, ; ] Aj > holds Vz. Also if x^ 2 , Xj-i, Xj, x i+1 are 
coplanar then Tj(t) = [nij_i Li mi\ = [L,-_i Li L i+ i] = 0. Thus coplanarity condition is also 
satisfied. 

14.2 Convexity and inflection preservation 

It satisfies convexity and inflection preserving criteria with suitably chosen tangent length. This 
is mainly due to the reason that apart from the tangent vector at a data point x, being coplanar 
with the plane IT the two consecutive sides Li and L i+ i lies on one side of it. 

14.3 Collinearity preservation 

Also if the data points Xj_i, x, ; and x i+1 are collinear and Lj_i • > then the tangent vector 
of the spline at x, ; is collinear with Xj_i, Xj and x i+1 . Thus the Catmull-Rom splines having 
tangent vectors with our magnitudes satisfies collinearity preservation criteria. 

We also observe that for collinear data arc, the shape of Catmull-Rom splines may not be 
aesthetically pleasing. We need to deviate from Catmull-Rom splines in accrodance with the 
conditions of modified definition of collinearity preservation criteria stated in sectional 

51 



Figure 22: By suitable choice of magnitudes of tangent vectors of Catmull-Rom splines can 
preserve the shape of data polygon 

15 Conclusion 

We have analyzed the characterization for shape preservation criteria for splines. We have 
improved upon the definitive criteria for convexity preservation for splines. We have studied in 
detail the inflection criteria and various results concerning it. We have also stated the results 
from the literature which in conjunction with our analysis are observed to give negative results 
regarding convexity and inflection preservation criteria for splines. Such negative results would 
have been difficult to perceive intuitively. We have also discussed the analysis for collinearity, 
torsion and coplanarity preservation criteria. 

We obtained a very important theorem 112.11 which states that there is a possibility that 
torsion preserving spline may need to have torsion set to zero at some nodepoints in order to be 
torsion continuous spline curve. From the literature we find that 

• Shape preservation criteria also gives better way of segmentation of curves. 

• Such curves can be very usefull tool for data reduction, which is very important for data 
transmission. 

• Such interpolation can be used for robot path determination. 
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We are currently working towards getting profound theoretical and experimental results in the 
above directions. 

We have found the characterization of all the shape preservation criteria for splines for the 
cubic case in terms of data points and slope vectors on them. 
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